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Abstract 

The de Sitter invariant Special Relativity (dS-SR) is a SR with constant curvature, 
and a natural extension of usual Einstein SR (E-SR). In this paper, we solved the 



£^ . dS-SR Dirac equation of Hydrogen by means of the adiabatic approach and the quasi- 

r~| ■ stationary perturbation calculations of QM. Hydrogen atoms are located on the light 

cone of the Universe. FRW metric and ACDM cosmological model are used to dis- 
cuss this issue. To the atom, effects of de Sitter space-time geometry described by 
Beltrami metric are taken into account. The dS-SR Dirac equation turns out to be a 
time dependent quantum Hamiltonian system. We revealed that: l,The fundamental 
physics constants me, e variate adiabatically along with cosmologic time in dS-SR 
QM framework. But the fine-structure constant a = e^jihc) keeps to be invariant; 
2,(2s^/2 _ 2pV2).spiitting due to dS-SR QM effects: By means of perturbation theory, 
that splitting A£'(z) were calculated analytically, which belongs to C'(l/i?^)-physics of 
dS-SR QM. Numerically, we found that when \R\ ~ {lO^G/y, Y^^Gly, W^Gly }, and 
z ~ {1, or 2}, the AE{z) » l(Lamb shift). This indicate that for these cases the 
hyperfine structure effects due to QED could be ignored, and the dS-SR fine structure 
effects are dominant. This effect could be used to determine the universal constant R 
in dS-SR, and be thought as a new physics beyond E-SR. 
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1 Introduction 



Einstein's Special Relativity (E-SR) has global Poincare- Minkowski space-time symmetry. 
E-SR indicates the space-time metric is r]^^ = diag{+, — , — , — }. The most general transfor- 
mation to preserve metric rj^j^i, is Poincare group. It is well known that the Poincare group 
is the limit of the de Sitter group with the sphere radius i? — j- oo. Thus people could pursue 
whether there exists another type of de Sitter transformation with R — finite which also 
leads to a Special Relativity theory (SR). In 1970's, Lu, Zou and Guo suggested the Spe- 
cial Relativity theory with de Sitter space-time symmetry (dS-SR) [I] [2]- In recent years, 
there are various studies of this theory [3]. In 2005, Yan, Xiao, Huang and Li performed 
Lagrangian-Hamiltonian formulism for dS-SR dynamics with two universal constants c and 
R, and suggested the quantum mechanics of dS-SR [1]. There is one universal parameter 
c (speed of light) in the Einstein's Special Relativity (E-SR). By contrast, there are two 
universal parameters in the de Sitter Special Relativity (dS-SR): c and R (the radius of de 
Sitter sphere and to character the cosmic radius). In this present paper, we try to study one- 
electron atoms, typically Hydrogen atom, of a distant galaxy (e.g., a Quasi-Stellar Object 
(QSO)) by means of dS-SR Quantum Mechanics (QM) suggested in Ref. [1]. 

As is well known that one of GR principles is existence of Locally Inertial System (LIS) at 
any point with small enough vicinity region in the curved space-time. In LIS, the expressions 
of physics laws are the same as ones in SR. Therefore determining the energy level shifts of 
a distant Hydrogen atom due to dS-SR QM will be useful to the cosmology when the curved 
space-time is the Friedmann-Robertson- Walker (FRW) Universe. 

Ref. [1] shows that the dS-SR dynamical action for free particle associates the dynamics 
with time- and coordinates-dependent Hamiltonian. In other hand, the Noether theorem 
assures the symmetry's Neother charges to be conserved even though that the Hamiltonian is 
time- and space-dependent. In Ref. [1], 10 external conserved Neother charges for dS-SR have 
been explicitly presented, which are free particle's energy, 3 momenta, 3 angular-momenta 
and 3 boost generators (see, Eqs (52)- (56) in [1]). Thus, the energy conservation law in 
dS-SR holds, and at the same time the dS-SR dynamics is a time-dependent Hamiltonian 
system. Contrasting with E-SR dynamics, this is a remarkable feature of dS-SR. This will 
cause time-dependent level shifts in atomic physics, and lead to some remarkable observable 
effects in cosmology, for instance, the physics constants varying adiabatically [9] and some 
specific level shifts for Hydrogen atom caused by time interval bing on cosmic scale. In this 
paper, we will focus on the splitting effect between 2s^^'^ and 2p^/^ states of Hydrogen in 
dS-SR Quantum Mechanics (QM). 

In this paper, the adiabatic approach [5] [6] [7] [8] will be used to deal with the time- 
dependent Hamiltonian problems in dS-SR QM. Generally, to a H{x,t), we may express 
it as H{x,t) = Hq{x) + H'{x,t). Suppose two eigenstates \s) and \m) of Hq[x) do not 
generate, i.e., = hium — ojs) = fu^ms 7^ 0. The validness of for adiabatic approximation 
relies on the fact that the variation of the potential H'{x,t) in the the Bohr time-period 
{ATms°^^^)H'^g = [211 / ujms)H'^g is much less than fiWras, where H'^^ = {m\H'{x,t)\s). That 
makes the quantum transition from state \s) to state |m) almost impossible. Thus, the non- 
adiabatic effect corrections are small enough (or tiny) , and the adiabatic approximations 
are legitimate . To the wave equation of dS-SR QM of atoms discussed in this paper, 
we show that the perturbation Hamiltonian described the time evolutions of the system 
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H'{x,t) oc {c^t'^/R'^) (where t is the cosmic time). Since R is cosmologically large and 
R » ct, the factor [c^t^ / R^) will make the time-evolution of the system is so slow that 
the adiabatic approximation works. We shall provide a calculations to confirm this point 
in the paper. By means of this approach, we solve the stationary dS-SR Dirac equation for 
one electron atom, and the spectra of the corresponding Hamiltonian with time-parameter 
are obtained. Consequently, we find out that the electron mass me, the electric charge e 
and the Planck constant h vary as cosmic time going by, but the fine structure constant 
a = e^Hhc) keeps to be invariant. Those are interesting consequences since they indicate 
that the time- variations of fundamental physics constants are due to well known quantum 
evolutions of time-dependent quantum mechanics that has been widely discussed for a long 
history (e.g., see [7] and the references within). 

The life time of a stable atom, e.g., the Hydrogen atom, is almost infinitely long. We 
can practically compare the spectra of atoms at nowadays laboratories to ones emitted 
(or absorbed) from the atoms of a distant galaxy, typically a Quasi-Stellar Object (QSO). 
The time interval could be on the cosmic scales. Such observation of spectra of distant 
astrophysical objects may encode some cosmologic information in the atomic energy levels at 
the position and time of emission. As is well known that the solutions of E-SR Dirac equation 
of atom are cosmologic effects free because the Hamiltonian of E-SR is time-independent, 
and the solutions at any time are of the same. Thus, after deducting Hubble red shifts, any 
deviation of cosmology atom spectrum observations from the results of E-SR Dirac equation 
could attribute to some new physics beyond E-SR. dS-SR is one of the most straightforward 
answers to such kind of deviations. 

In E-SR Dirac equation of Hydrogen, 2s^/^- and 2p^/^-states are completely degenerate. 
The hyperfine effects of QED break this degeneracy, and turn out famous "Lamb Shift": i.e., 
^2,1/2 - E2pi/2 = 1 Lamh Shift = 4.35152 x IQ-^eV = 1057.9MHz x 2'Kh. In this paper, 
we solve the dS-SR Dirac equation for one-electron atom and reveal a dS-SR effect which 
also contribute a level shift to break the 2s^/^ — 2p^/^-degeneracy. This shift is proportional 
to Q^Q^/R^, where is the distance between the Earth and an observed galaxy (e.g., a 
QSO), is the corresponding time interval, and R is the radius of de Sitter sphere. When 
Q^Q^ /R^ were large enough, the 2s^/^ — 2p^/^ splitting due to dS-SR effects would be much 
larger than QED's Lamb shift, the observation of this splitting in cosmo logical experiments 
could provide a criteria to check dS-SR. 

The contents of the paper are organized as follows: In section II, we briefly recall the 
classical mechanics of de Sitter special relativity and the corresponding quantum mechanics. 
The dS-SR Dirac equation for spin-1/2 is presented; Section HI is devoted to discuss Hydro- 
gen atom described by dS-QM and embedden in light cone of Friedmann-Robertson- Walker 
(FRW) Universe; Section IV shows the solutions of usual E-SR Dirac equation for Hydrogen 
atom located at distant galaxy. Especially, the wave functions and energy values of the 
states 2s and 2p^/^ are presented explicitly. The Hydrogen atom energy level shifts due to 
gravity in FRW Universe are estimated; In Section V, we derive dS-SR Dirac equation for 
Hydrogen atom up to terms being proportional to l/R^] Section VI: dS-SR Dirac equation 
for spectra of Hydrogen atom; Section VII: Adiabatic approximation solution to dS-SR Dirac 
spectra equation and time variation of physical constants; Section VIII: 2s^/^ — 2j»^/^ split- 
ting in the dS-SR Dirac equation of Hydrogen; Finally, in Section IX, we briefly discuss and 
summarize the results of this paper. In Appendix A, we derive the electric Coulomb Law in 
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QSO-Light-Cone Space; In Appendix B, we show the calculations of adiabatic approximative 
wave functions in dS-SR-Dirac equation of Hydrogen in detail. In Appendix C, we provide 
analytic calculations to the matrix elements of perturbation Hamiltonian in dS-SR, which 
yield (25"^/^ — 2p^/^)-hyperfine splitting discussed in the text. 



2 Special Relativity with de Sitter Symmetry and dS- 
SR Dirac equation 

The precise dS-SR theory were formulated in 1970-1974 by LU, ZOU and GUO [I] [2] 
(for the English version, see, e.g., Ref. |1] [3]). Two theorems were proved: 

Lemma I: Inertial motion law for free particles holds to be true in the space-time char- 
acterized by Beltrami metric 

B^^{x) = + ^, .^ Vf^xVupX^x^, (1) 

where a{x) = \ — -^rj^^x^x'^ , and the constant R is the radius of the pseudo-sphere in de 
Sitter (dS) space, i?^ > or < that corresponds to dS symmetries 5'0(4, 1) or S'0(3,2) 
respectively. This claim means that in the space-time characterized by B^^, the velocity of 
free particle is constant, i.e., 

= if^ = constant, for free particle (2) 

which is exactly the counterpart of E-SR's inertial law in Minkowski space characterized by 
rjfj^i^. (see Refs. [3] |1] for the English version of proof to Eq.([2])). 

Lemma II: The de Sitter space-time transformation preserving B^y{x) is as follows 

x^'^i^' = ±a(a)^/V(a,a;)"^(x"-a")D^, (3) 
= L'^ + R~\ya\aia) + a^/^ia))'^L''^, 
L : = (L;;)g 50(1,3), 

(t(x) = 1 - —7]^yX^'x'', a{a,x) = 1 - —r]f,^a''x'' . 

where x'^ is the coordinate in an initial Beltrami frame, and x^ is in another Beltrami 
frame whose origin is in the original one. There are 10 parameters in the transformations 
between them. Under the transformation Q, we have the equation preserving 5^,^ as follows 

B,Ax) ^ B,^{x) = ^^B,,ix) = B,^{x). (4) 

( see Appendix of Ref. [1] for the English version of proof to Eq.(j4j)). Eq.(jl]) will yield 
conservation laws for the energy, momenta, angular momenta and boost chargers of particles 
in dS-SR mechanics Here, we like to address that the space-time symmetry described by 
Eqs. (IH]) (jlj) is a global symmetry since both and are constants instead of functions of 
space-time x^. This situation is same as E-SR's, where the Poincare symmetry is global. 

Based on the dS-SR space time theory described in above two lemmas, the dS-SR dy- 
namics described by the Lagrangian as follows |1] 



ds ^/B^x)dx^ /n f \-a-u rr\ 

-moc— = -moc^ — = -mQC\l B^y[x)xf'x'' , (5) 
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where = -^x^, B^y[x) is Beltrami metric ([T]). Then the canonic momenta and canonic 
energy (i.e., Hamiltonian) reads 



TT,- 



H 



dL_ 

dx^ 



-m.oa(x)TBinX^ 



dL 

^ di^^' - ^ = moCcr{x)TBof,x'', 



(6) 
(7) 



i=l 



where 



ds 



1 



rjijX^x^ 



) + 2triijX^x^ — rjiji^xH"^ + 



{riijX'-x^] 



From the equation of motion = 0, we have [1] 



|x'=0 



0, 



(9) 



whose corresponding one in E-SR is 



7\xi=0 



d 



1 — v"^ /c^ 



0. 



(10) 



t)=constant 



By means of the standard procedure to perform the canonic quantization, we obtained the 
dS-SR wave equation for spinless particle 



2^2 



0, 



which is just the Klein-Gordon equation in curved space-time with Beltrami metric B^y. 
The measurable conserved 4-momentum operator is [1] 



ih 



5x^ 
2^ 



(12) 



The corresponding Dirac equation which describes the particle with spin 1/2 [10] [12] reads 



(13) 



where is the tetrad and is the covariant derivative with Lorentz spin connection u'^. 
Their definitions and relations are follows (e.g., see [12]) 



d^--co:''aab, D^ = B^^D, 



ab 



1 1 

2?7"^ OTab = -[la, lb] , - Wab, (^cd] = Vacf^bd " Vadf^bc + VbdC^ac " VbcCTad, 



Bf'idxBy. + d.ByX-dyB, 



(14) 



It is straightforward to check that the components ^/'q, (a = 1, ■ ■ ■ 4) of the spinor satisfy the 
Klein-Gordon equation ffTTl) . 
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3 Hydrogen atom embedded in light cone of Friedmann- 
Robertson- Walker Universe 



The isotropic and homogeneous cosmology solution of Einstein equation in GR (General 
Relativity) is Friedmann-Robertson- Walker (FRW) metric. In this section we discuss the 
Hydrogen atom embedded in FRW Universe and described by dS-SR Dirac equation. 

As is well known that GR can be viewed as a kind of gauge field theory [10] [11]. The 
dynamic equation of GR can be yielded by means of the localization of external global space- 
time symmetries, such as Lorentz group SL{2,C), transition group T4, Pioncare group etc. 
Such localizations make global space -time transformation — )■ x''^ = K^^^x^ + to be x^ — >■ 
x'M = A.^^^{^x)x'^ + a^{x) = f^{x) which is arbitrary curvilinear coordinate transformation with 
Christoffel symbol connection {^^} for the torsion-free space-time, and leads to construct GR 
by means of Riemann tensor. It is essential that the the framework of GR is generally free 
to the gauge theory's underline global external space time symmetry for the torsion-free 
space-time. Now let's see the global de Sitter space-time transformation Eq.Q. When 

— )■ a'^(x), — )■ the transformation is localizated to be x^ ^ x^ = f{x)^ which is 

an arbitrary curvilinear coordinate transformation. The connection is still Christoffel symbol 
and no corrections are yielded to the GR framework, and hence the considerations of atoms 
described by dS-SR QM in FRW Universe are legitimate. 

One way to detect the spectrum of distance atom is spectroscopic observations of gas 
clouds seen in absorption against background Quasi-Stellar Objects (QSO), which can be 
used to search for level shifts of atom for various purposes (see, e.g., [13] [II])- In the 
observations of gas-QSO systems in the expanding Universe, one can observe two species of 
frequency changes in atomic spectra: the Hubble redshift (z) caused by the usual Doppler 
effects and a rest frequency change due to the dynamics of atom beyond E-SR. The latter can 
be found by measuring the relative size of relativistic corrections to the transition frequencies 
of atoms on the gas-QSO (or on QSO for briefness). A widely accepted assumption is that 
this rest frequency change is independent of the Hubble velocity and cosmologic acceleration 
of the gas in the Universe. Therefore all relativistic correction calculations for atom spectra 
were performed in a "rest" inertial reference frame without any Lorentz boost and non- 
inertial effects caused by frame-origin motion [T5|[T6] . In this present paper, the calculations 
based on dS-SR Dirac equation are performed in such rest reference frame. In this framework, 
the time-varying (2s^/2 - 2p^/2)_gpijt^jng /\e{2s^/^ - 2p^/'^) = AE{t) will be calculated. 
Furthermore, the t — z relation has been established for the Universe with the observed 
acceleration in ACDM model [18l[20]. Employing this relation, we have AE{t) = AE[t{z)). 

Now, we show the dS-SR Dirac equation of Hydrogen atom on a QSO in the Earth-QSO 
reference frame. As illustrated in FIG.l, the Earth is located at the origin of frame, the 
proton (nucleus of Hydrogen atom) is located at Q = = ct, Q^, = 0, = 0}. To 
an observable atom in four- dimensional space-time, the proton has to be located at QSO- 
light-cone with cosmic metric (7^^. Namely, Q must satisfy following light-cone equation (see 
FIG. 1 and set = = for simplification) 

ds^ = g^,u{Q)dQ^dQ'' = 0, (15) 

which determines = f{Q^)- We emphasize that the underlying space-time symmetry for 
the atom near Q described by dS-SR dynamics is de Sitter group instead of to limit it as 
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Poincare symmetry of E-SR as usual, which is only a special limit of dS-SR's. The corre- 
sponding space-time metric is B^,^{Q) = r]^^ ^1 + - W ) j _^ _^^^^gA^^^gp _^ 0{1/R^) 
(see Eq.([T])). Note B^y{Q) is position Q-dependent, and Lorentz metric 1]^^ for E-SR is not. 
Explicitly, from Eq.(IT]), we have 



/ . 2(Q0)2-(Q1)^ 

' ^ 7?2 



V 

' ^ 7?2 



-1 + 



ij2 



-1 - 




\ 






2(Q^)^-(Q°)^ 






-1 - 



-1 





(Q0)2-(Q1)2 






(Q0)2-(Ql)2 
R2 





-1 + 






(Q0)2-(Ql)2 , 
i?2 / 

o\ 




(Q0)2-(Q1)2 



(16) 



(17) 



i?2 



In following, we will solve Eq.( fT5ll to determine in FRW Universe model. 



Qt I Earth 




Electron 



Nucleus 



Figure 1: Sketch of the Earth-QSO reference frame. The Earth is located in the origin. The position 
vector for nucleus of atom on QSO is Q, and for electron is L. The distance between nucleus and 
electron is r. 



The Friedmann- Robertson- Walker (FRW) metric is (see, e.g., [T7] 



c^dr - a{ty 



1 — kr'^ 



r^d6^ + sin^ 



= {dqy -a{ty {dQ'dQ' + 

= 9>..{Q)dQ^dQ' 



KQ'dQ'f 
1 - kQ'Q' 



;i8) 
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where r = \/Q^^, = r sin 6' cos 0, = r sin 6' sin 0, = rcosO has been used. As is 
well know FRW metric satisfies homogeneity and isotropy principle of present day cosmology. 
When Q'^ = = 0, from f|T8l) . we have 

/ 1 HQ'^Y 
9iiu{Q) = V/^iy - a{tf6f,i6^i -—r + 1 + 



a(t)2 l-k{Q^] 
-ia{t)^-l)i6^2Su2 + SM- (19) 

For simpleness, we take k = and a{t) = 1/(1 + z{t)) (i.e., a(to) = !)• And the red shift 
function z{t) is determined by ACDM model [l8H20](see, e.g., Eq.(64) of [19]): 



dz' 



H{z'){l + z'y 



(20) 



where 



Ho = 100 /i ~ 100 X 0.705km ■ s'^/Mpc, 
n^o ^ 0.274, nnor^lO-'. (21) 



Figure of t{z) of Eq.([20]) is shown in FIG.2. 




Figure 2: The t - z relation in ACDM model feQ. ((20)) ). 

From (fT9|) . the FRW metric reads 

Qt^uiQ) = VtMu - {a{ty - 1) (5^1 5^1 + 5f,2S^2 + 5^35^3)- (22) 
Substituting fl22l) into f|T5|) . we have 



rfgo = J^rfgi = a(t)rfgi= \. dg\ (23) 



Consequently, by using Eq.f l20|) and = ct, we get desirous result: 



e'-Tj^. (24) 
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Figure of Q^{z) of Eg. (12^ is shown in FIG. 3. Ratio of over Q° is shown in FIG. 4. 

Then the location of distance proton is {Q^, Q^, 0, 0} in the space-time with FRW metric. 

We treat Hydrogen atom as a proton-electron bound state described by quantum me- 
chanics under instantaneous approximations (see FIG. 1). The electron's coordinates are 
L = {L° = ct^ — ct, L^, L^, L^}, and the relative space coordinates between proton and 
electron are = L* — Q*. The magnitude of r = ^J —riijX^x^ ~ as (where — 0.5 x 10~^°m 
is Bohr radius), and ~ a^. 

According to gauge principle, the electrodynamic interaction between the nucleus and 
the electron can be taken into account by replacing the operator in eq. (fT3|) with the 
f/(l)-gauge covariant derivative = — ie/{ch)A'^, where A^^ = {0Bj A}. Hence, the 
dS-SR Dirac equation for electron in Hydrogen at QSO reads 

{te,a7^V^ - = 0, (25) 

where /i = me/(l + ^) is the reduced mass of electron, = ^| jco^^'^'-aab — ie/{ch)A'^, 

and have been given in eqs.(IT]) ( IT^ . For our purpose, we approximately write and 
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uf up to as follows: 

e^a = {^-^^)^'a-'^^ + 0{l/R% (26) 

= ^KL'-vtn + Oil/R'). (27) 

In the following sections we are going to solve dS-SR-Dirac equation for Hydrogen atom 
on QSO in FRW Universe model. In this quantum system, there are two cosmologic length 
scales: cosmic radius, say R ~ lO^^/y, and the distance between QSO and the Earth, that 
is about ~ ct: say R > ct > lO^ly, and two microcosmic length scales: the Compton 
wave length of electron Qc = h/{mec) ^ 0.3 x 10~^^m, and Bohr radius = ^^/(meC^) ~ 
0.5 X 10~^°m. The calculations for our purpose will be accurate up to O^cH"^ / R^). The 
terms proportional to 0{cH'^/R'^), 0{ctac/R^), 0{ctaB/R^) etc will be omitted. 

4 Solution of usual E-SR Dirac equation for hydrogen 
atom at QSO 

4.1 Eigen values and eigen states 

At first, we show the solution of usual E-SR Dirac equation in the Earth-QSO reference 
frame of Fig.l, which serves as leading order of solution for the dS-SR Dirac equation with 
i? — oo in that reference frame. For the Hydrogen, — )■ = — ie/{ch)A'^j (noting 
uj"'^'^\r^oo = 0), where A^.j = {(pM^x), Am{x)}, and (f)M{x) and Am{x) are nucleus electric 
potential and vector potential at in Minkowski space defined by following equations (see 
Appendix A, Eq. ffT39|) ) 

-7]'^d,d,Mx) = V20m(x) = -47rp(x) = -47re(5(=^)(x), (28) 

V{dxAlj)-d^AM = -— J = 0. (29) 

The solutions are (Pm{x) = e/r and Am = 0. And hence do — j- Vq = do — irjooe^ / (chr) . 
Then, the E-SR Dirac equation reads 

ihdtip = (^-ihca ■ Vl + /ic^/? - (30) 

where /3 = 7°, a* = PY- Noting the nucleus position Q =constant, we have 

= I: = diiw) =¥r^^' ^^^^ 

and eq. (l30l) becomes the standard Dirac equation for electron in Hydrogen at its nucleus 
reference frame. Energy W for E-SR-mechanics is conserved (hereafter, we use notations 
of [22]), and the Hydrogen is the stationary states of E-SR Dirac equation. The stationary 
state condition is 

ihdtip = Wi/j. (32) 
As is well known, combining eqs. flSU]) . f lHT]) with f l5^ . we have 

Wi/j = (-ihcd ■ V + /xc^/S -—']'ip = Ho{r, h, e)i), (33) 
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which is the stationary E-SR Dirac equation for Hydrogen. The problem has been solved in 
terms of standard way, and the results are follows (see, e.g., |21] |22j ) 

« = -, |«:| = (j + l/2) = l, 2, 3 ■■■ 
s = ^ — a^, 77, = 1, 2, 3 ■ ■ ■ . 
And its expansion equation in a is 

^-'"'^-^^^'A^-^^y-)- '''' 

The corresponding Hydrogen's wave functions i\} have also already been finely derived (see 
e.g., The complete set of commutative observables is {if, k, j^, j^}, so that = 

^n,K.,j,j^{r, K fJ', ce) = ?/;™-'(r), where j = L + |5", hn = f3{a ■ L + /l), and a = e^/{hc). The 
expression of ip^^ (r) is as follows 



where 



^rW=( ^f'i^f^ul ] (36) 



1/2 

X™^(f) = E Ci;:^-^^^,^.^Yr-'^^{r)r^ (37) 

ms=-l/2 



sm ( 



and 



g^{r) = 2Xikc + WcY/'e~''-{2Xry-'a', 

X (^nM{l - n, 2s + 1, 2Ar) + - M(-n', 2s + 1, 2Ar)^ , (39) 

/.(r) = 2X{kc-WcY^'e~'^{2Xrr~'a', 

X |^n'M(l - n', 2s + 1, 2Ar) - (^^ - M(-n', 2s + 1, 2Ar)^ , (40) 

where 

kc = fic/h, Wc = W/{ch), A = (4 - Wl,y^^, s = Vk^ - n' = n - \k\, 

n = -ijij + l)-lil + l) + l/A), (41) 

where W has been given in Eq.( !34l) . and M(a, b, z) is the confluent hypergeometric function: 

Ma,6,2; =l + x + yf^ + ■■■^7^ + ■■■ 
6 2!(6)2 n\{b}n 

(a)o = 1, (a)„ = a{a + l)(a + 2) • ■ ■ (a + n - 1), 
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and ol is the normalization constant required by 

/■oo 

Jo 

To Hydrogen-like one electron atoms with Z > 1, the energy level formula and the eigen 
wave functions expressions are all the same as Eqs. (j34l) - (l4Tl) except a ^ = Za. In FIG. 5 
the levels of one electron atom with Z = 92 are shown. 











n 1 


j *■ 


a 


W— Hi,,c (eV) 


"2 1 


3/2 -2 


4 


-29664 


-2 0,1 


1/2 -1,1 


4 


-34235 


J 


1/2 -1 


2 


-132353 











Figure 5: The energy levels of the one-electron atom with Z = 92. The figure shows relativistic energy levels 
calculated using equation p4|) with replacement of a — > ^ = Za. They are labelled by the quantum numbers 
n, I, j, K and their degeneracy d, up to n = 2. 



It is learned from above that since R — oo, Bf^„ — 77^^, dS-SR— )-E-SR, the Hamiltonian 
of E-SR is cosmic time independent. So, the spectra of Hydrogen at any place in FRW 
Universe are the same, and there is almost no cosmology information of the Universe in the 
spectrum solutions Eq. fl51|) of E-SR Dirac equation. 

4.2 2s^/^ and 2j9^/^ states of Hydrogen 

As is well known that the state of 2s^/^ and state of 2j9^/^ are complete degenerate 
to all order of a in the E-SR Dirac equation of Hydrogen. Namely, from (l34l) and k, = 
+ 1) - /(/ + 1) + 1/4), we have 

AW{2s'/' - 2p'/^) = W^n=2, .=-1) - W^(n=2, .=+1) = 0. (42) 

By means of Eqs. fl36l) - (l40l) . the wave functions of states 2s^/^ with k = —1 are as follows [22] 

.m, f(2.i/2) (r)x;r'(f) (2,1/2) \ , . 

\ ^/(2si/2)(,rjX-Kirj(2si/2) / 

where 

, ^ l{2Xy'+'kc{2s + l){kc + Wc) fWc A / , 2Wc\ \.., 

^^^--^^^^-^ 2Wci2Wc + kc)Ti2s^l) ^ ' (l^ " 2in + ^ j T'^ 

V 2Wai2Wc + kc)Ti2s + l) ' ' ' 2^1^^^)' 
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where the expressions of notations kc, Wc, A and s have been shown in Eq.( HT]) following 
[22], and due to Eqs.dST]) ([38]) 



x;^/'(f)(2.v.) = ( ^f), X."\i)i2s^i^) = f °o ) ' ^^th = (46) 



1/2,.. / -cos^Fo" \ / -sin^e-^^yo" \ 

For the 2p^/^ state with k' = 1 



(2P)J-1/2V V «/(2pi/2)(^)X_«'(r)(2pV2) ' 



where 



. ^ (2A)2-+^fcc(2. + l)(fcc + H^c) _,,/Wc . ^ A / 2iyc\ w... 



, , (2A)^-+^fcc(2. + l)(fcc - Wc) , UWc A / 2iyc\ \ 

^^^^^^^^^'^-1 2Wc{2Wc-kcn2s + l) ' ' ^^ + 27nl'"^JV 

-1/2/. N _ ( \fz^^ 



(50) 



X.. (f)(2pV2) = I I ' (^)(2P^/^) = I V^n,. I ' (51) 



X-K'(r")(2pi/2) 



i cos er^oj^^^) - y I sin ee~''t'Yl 
sin ee*<^y 0(^0) + ^co^eV^ 



^'^^^^^"^ " I ^sin^e^^n-(^0) + ,[\co.eYX ^ ' ' 



where Y^{9(t)) = Jf cos^, F/(^0) = -Jf e^-^sin^, Ff ^(^0) = e-^-^sin^ 



4.3 Hydrogen atom energy level shifts due to gravity in FRW 
Universe 

Secondly, we estimate the influence of external cosmological gravitational fields on the 
Hydrogen energy levels. 4D-FRW space-time is curved, and the atom locates in its tangent 
flat space-time described by B^^\ji^ao = V^iu■ It is well known that the interaction between an 
external gravitational field and atom may cause some energy level shifts of the atom [23] [2l] . 
In [21], gravitational perturbations of the Hydrogen atom are derived by means of Fermi 
normal coordinate method [25]. It is found that energy level shift of one electron atom to 
the first order of curvature for nS states (see Eq.(8) in 



E^'\nS) = Ic" V'n'(5n' + l)7^oo, (53) 

where fi~^(~^ = a^, C = and TZqo is Reimann curvature tensor. By means of FRW metric 
fl22|) . Einstein equation (e.g., see Eq. (1.5.18) in pp.36 of [IT]) and ACDM, we have 

7^oo = -3d/{c^a) = —Gpdrtmo - 2Qa) ~ -1.05 x 10~^^cm-\ (54) 
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and hence 

E^'^\nS) = --n\5n^ + 1) x 7.3 x 
6 

= -ln^{5n^ + 1) X 1.7 x 10^^\Lamb shift). (55) 
6 

It is an extremely tiny level shift. To other energy levels, the corresponding level shifts are 
about same order magnitudes of Eq. (!55|) . Physically, it is noting but a cosmologic gravity 
tide effects on the atom described by E-SR Dirac equation. Eq. (ISSi) indicates that such tide 
effects can be completely ignored. Namely the Coulomb potential in the atom is determined 
by Eqs.([28]) ([29]) which arise from Eq. f ll3Tl) of Appendix A with g^j^i, — "Ufiui and any influences 
of cosmological FRW metric are able to be ignored. Similar conclusion were reached in the 
de Sitter Universe model in [23]. In this present paper, though the atomic dynamics is of 
dS-SR Dirac equation rather than E-SR's, we will also ignore this sort of cosmological gravity 
tide effects to atoms. Namely, instead of g^^ = rj^y for E-SR QM, we shall employe Beltromi 
metric g^y = B^^ to derive the Coulomb potential in the atom for dS-SR QM in the follows. 



5 dS-SR Dirac equation for Hydrogen atom 

By eqs.OSS]), Oil]), (127]), 9^ ^ = ^ - {u'^^'^aab - ie/{ch)At', and noting = B^^Vl, 
we have the dS-SR Dirac equation for the electron in Hydrogen at the earth-QSO reference 
frame as follows 



2i?2 



^Mpi 



h 



^ = 0, 



(56) 



where factor hcP in the front of the equation is only for convenience. We expand each terms 
of (l56l) in order as follows: 

1. Since observed QSO must be located on the light cone, then VabL^'L^ ~ VabQ'^Q'' = 
{Q^y — {Q^y, and the first term of (156]) reads 



hcl3i 1 



2i?2 



2i?2 



with /37^ = {/37° = (3^ = 1, (3f = a'} 



hc(3i-f''Vl:^p (57) 
(58) 



hcPtYVj:^lj = {ihdt + hca-V + -^oofYcTab + e^Y'B^yA'')^/;, (59) 

where A'^ = {A^ = <Pbi A^} are determined by Maxwell equations within Constance 
metric g^y = B^y{Q) and f = {f = cpj^^^, f = 0} (see Appendix A): 



-B'^{Q)did,Mx) 



1 - 



E? 
—Awe 



1\2 



{Q 



1\2 



i?2 d{x^)\ 



^b{x) = J 

c 



v/-det(5,,(g))5oo(g) 



5(3)(, 



Air 



d'd.A'' - B'^'d.dyA' = f = 0. 



(60) 
(61) 
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The solution is (see Appendix A) 



3(g°)^-2(g 



(62) 



where tb = \/ (5^)^ + (x^)^ + [x^Y with = + ^^s-]- In the follows, we 

use variables {a:^,x^,x^} to replace {x^, x^, x^}. Following notations are introduced 
hereafter: 



d 



d 



d 



x' G {x^, x^, x^} 



Sx-"^ dx"^ dx^ ' 
Then, noting ^ = fi J, = J, - J,, the eq^ becomes 



a 



2i?2 dx^ 



iAio't + ihca ■ \ B — iiic a 7-— + 



,a,b fi 



+ 



i?2 



1n2 



2i?^ " 9x1 ' 4 7 



(63) 
(64) 



(65) 



2. Estimation of the contributions of the fourth term in RSH of ( l65l) ( the spin-connection 
contributions) is as follows: By f l27|) . the ratio of the fourth term to the first term of 
am is: 



he ct 



ct h 



1 Ctttr 



4 2i?2 m,c2 8i?2 8 i?^ 



0, 



(66) 



where Cc = h/ (mgc) ^ 0.3x lO^^^m is the Compton wave length of electron. 0{ctac/ B?)- 
term is neglectable. Therefore the 3-rd term in RSH of (159|) has no contribution to our 
approximation calculations. 

3. Substituting (ES]) (ES]) ([58]) into ([57D and noting L"" ~ Q'^ (see FIG. 1), we get the first 
term in LHS of (!56l) 



1^2N / Q 

ih— + ihca ■ Vs 



dt 



-mc-rrrT^a -— + [1 H ^^^5 Je0B ^^" 605 ) ■j/' 



2/22 dx 
o\2 ,'/nl^2 



i?2 



/22 



2/?2 



+ 



1 + 



2i?2 



1\2 



e(pB 



a^e(f)B + 0{l/R^) }iJ.{67) 
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4. The second term of is 



I 

2^2- 



ihc 



26 



en 



(68) 



where following approximation estimations were used 

(L') {L') as „ 
R R R 

Noting L° ~ Q°, ~ Q\ (EH]) becomes 



2i?2 



2i?2 



a ihdf 



{Q 



0\2 



2/22 



-the—— H ——tnca 



2i?2 2i?2 



(69) 



2i?2 



(70) 



5. Therefore, substituting (1671) ( ITOl) into ( l56l) . we have 



V 2i?2 2/^2 y 



—ihc 



U^2 



or 



2i?2 



2i?2 



(Q 



0^2X g2 



2i?2 / 



w —a 

^ 2i?2 



+ 



2i?2 dx^ 



(71) 



(72) 



3(Q°)^-2(Qi) 



2^^^^)e/rB were used (see Eq.([62])). Eq.([72]) is dS-SR Dirac 



where 0b = (■ 

wave equation to the first order of 0{j^). Two remarks on (1721) are as follows: 

i) When R — )■ oo, Eq. (172|) goes back to usual E-SR Dirac equation of Hydrogen, which 
has been discussed in the last section. 

ii) Eq. (l72l) is a time-dependent wave equation. It is somewhat difficult to deal with 
the time-dependent problems in quantum mechanics. Generally, there are two approx- 
imative approaches to discuss two extreme cases respectively: (i) The modification in 
states obtained by the wave equation depends critically on the time T during which the 
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modification of tlie system's "Hamiltonian" take place. For tliis case, one would use 
the sudden approach; And, (ii), for case that of a very slow modification of Hamilto- 
nian, the adiabatic approach works [6]. To wave equation of fl72|) . like the discussions 
in Introduction of this paper, since \R\ is cosmologically large and \R\ » ct, fac- 
tor {((5°)^/-R^, {Q^)'^/R'^} oc [cH'^/R^) makes the time-evolution of the system is so 
slow that the adiabatic approximation |5] may legitimately works. In the below (the 
subsection E), we will provide a calculations to confirm this point. 



6 dS-SR spectra equation of hydrogen 

In order to discuss the spectra of Hydrogen by dS-SR Dirac equation, we need to find out 
its solutions with certain physics energy E. By eq.f ll2p . and being similar to f l32|) . the dS- 
SR-energy eigen-state condition for fl72]) can be derived by means of the operator expression 
of momentum in dS-SR f lT^ : 



E 
p = — 
c 



E 

Eijj 



ih 
ih 
ih 



lyf^L , Set 



dt- 
1 - 



dt 



hct 



0\2 



i?2 



(73) 



where a estimation for the ratio of the 3-rd term to the 2-nd of Eijj were used: 



"2/12 • 



\ih 



2R? 



5h 



i?2 



ihdtip\ 



-2c2t2 
ij2 



^1 



5 Oc 
2ct 



where Oc — 0.3 x 10 ^^m is the Compton wave length of electron and ct is about the 
distance between earth and QSO. In our approximative calculations ac/{ct) is neglectable. 
For instance, to a QSO with ct ~ lO^ly, ac/{ct) ~ 10"^^ << {ctf/R^ ~ 10^^ Hence the 
3-rd term of Eip were ignored. 

Inserting fl73|) into (172|) . we obtain the 57?.c_R-spectra equation of hydrogen 



1 + 



0\2 



Eip 



a\2 



2i?2 

2i?2 



ihca - V B + I 1 



my - ( Q 

2i?2 



1\2 



fic'^/3 



e 

tb 



-ihcT^ -Vb + /ic'^/S 

rB 



■ihc- 



2i?2 
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or 



0\2 



- I 1 

+ 



ihcd ■ Vb + ( 1 



1\2 



2i?2 J rB 



2/^2 



-a 



—ihc^ ■ V 



the 



2i?2 dx^ 



+ 



-ihtcd -Vb + jJ'tC P 

Tb, 

inc- 



2i?2 



-a 



-ihcl^ -Vb + 

Tb 



2i?2 

= ( HoirB, ht, /it, et) + H' )tp = H^ds-SR)'ip, (74) 

which is up to 0{c^t^ /R^) (say again, C(1/_R^), 0{ctaB/ R^), 0{ctac/R^) terms have been 
ignored), and where 



HoirB, ht, fit, et) = -ihtcd ■ Vb + yUjc (3 



2n 



rB 



2i?2 



H' 



with 



et 



d \ 

-a^Ho{rB, ht, fit, et) + ^^c— j 

(Q0^2^ 



h-- 



and 



2i?2 
2i?2 
4i?2 



htc he 



4i?2 J 



e, 



a. 



(75) 
(76) 

(77) 
(78) 
(79) 

(80) 



Using notation in [22], Eq = W and the unperturbed eigen-state equation is 

Wijj = Ho{ht, /it, et)^ = (^-ihtca ■ Vb + /itc^/? - ^, (81) 

which is same as f l3^ except ^, /i, e being replaced by ht, fit, et- However, since the time t 
is dynamic variable in the time-dependent Hamiltonian system, at this present stage we do 
not know whether t can be approximately treated as a parameter in the system. Hence, we 
cannot yet conclude h, fi, e are time variations described by (!77|) . (fTHI) . (179|) at this stage. 
In the following section, we pursue this subject. 



7 Adiabatic approximation solution to dS-SR-Dirac spec- 
tra equation and time variation of physical constants 

Comparing (ISTl) with (l33l) . we can see that there are three correction terms in ( IHTl) . which 
are proportional to {cH'^/R'^) and service of effects of dS-SR QM. Since R » ct, we argue 
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that the corrections due to the effects should be small, and the adiabatic approach works 
quite well for solving this QM problem. In order to being sure on this point, we examine 
the corrections beyond adiabatic approximations in below by calculating them explicitly 
for a certain z. Suppose ^ ^ 3 ~ 4, FIG.(!4]) indicates ~ 1.7(5° = 1-7 ct, and hence 
(Q^)^ ~ 3((5°)^ = 3 c^t^ in fIST]) . Rewriting this spectra equation fIST]) in version of wave 
equation like eq.( l30|) via E =^ ihdt, we have 



ihdtiJ = H{t)tlj = [HoivB, h, /i, e) + H'^{t)]ilj, 



where 



Ho{rB,h,^,e) 



-ihcd - V B + l^c (3 (see eg. (133 



2R^ 



Ho{rB,h, 3/i,e). 



(82) 
(83) 

(84) 



Suppose initial state of the atom is ip(t = 0) = ipsivB.h, n.a) where s = {n^, k^, j^, jsz}, 
by eqs. (182|) (183!) (l84l) . and catching the time-evolution effects, we have (see Chapter XVII 
of Vol II of [6] , and Appendix B) 



l) ^m(rB,/it,/it,ef)e( ' 



Jo 



^5) 



where ipsi^B,fit, fJ't,^t) is the adiabatic wave function, ht, Ht are given in (177|) — (1791) . and 

-cH 



m\ {Ho{rB,h, 3yU,e) - Ho{rB, h, fi, e)) 

f^mi l^mi jm5 jmz\'^f^C (3\ng, Kg, jg, jsz)(^ ^ ^ 

m|/3|s)e-^('"^-'^'")*. 



-cH 



h 



^7) 



Note, formula {m\HQ{r,e)\s)\m^s = has been used in the calculations of ( l86l) . The second 
term of Right-Hand-Side (RHS) of eq. (ISSjl represents the quantum transition amplitudes 
from ips-staXe to ipm, which belong to the correction effects beyond adiabatic approximations. 
Now for showing the order of magnitude of such corrections, we estimate \HQ{t)ms/^ms\ 
for s = 1 = (Is^/^K = -l,mj = 1/2), m = 2 = (2s^/^ k = -l.nij = 1/2). Noting 
Wn ~ fic^ — fic^a"^ / (2n^), and the Compton wave length of electron Oc = h/ {niec) ~ h/ (fic) ~ 
0.3 X 10~^^m, from Eqs.(IH6])(IH7D we have 
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128 



(2|/3|1) 



with /3 



where the state (2| has been given in Eqs. 



- (ISUl) and the state |1) is as follows: 

1/2, 



^(isi/2)(r)x« (f)(i,i/2) 
2)(r)xiC(f)(i../.) 



(89) 
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where 



^(1.1/2)1'^) 



2kcT{2s + 1) 



'(2Ai)2^+i(^^ 



c 



W, 



2kcT{2s + 1) 

where Wci — meC^(l - a^/2)/{ch), Ai = a/A;^ - W^^ and 



(90) 
(91) 



Xy^(f)(l,l/2) 







, X-K (f)(lsl/2) 



cos er^ 



sm 



(92) 



By means of expressions of ^j'^ll'^l'^^ir) (Eqs.^-^) and V'[^I)"^ii/2(r) (Eqs.(l43D-(l50D), we 
have 



(2|^|1)= / cirr^ (^(1,1/2) (r)^(2,i/2)(r) - /(i,i/2)(r)/(2,i/2)(r)) 
Substituting ( 193|1 into ( IHSjl . we obtain 



-1.12 X 10' 



(93) 



5.6 X 10^ X 



ttc Ct 

RR 



(94) 



Considering Compton wave length of electron = — 0.3 x 10 ^^m ^ 0.3 x 10 "^^ly and 
both Q° = ct and i? are cosmological large length scales, and R > Q^, therefore we have 
that 



^21 



'1.7 X 10^24; go 
^ — X — << 1. 

R R 



(95) 



To generic (m| and \s), similar to the calculations of Eg. (1941) . we always have 

-^o(^)ni 



(constant) ■ ^ x 

R R 



Since the {constant) at last is about ~ 10^°, and then 
obtain 

'Kit), 



R 



R 



(96) 



we finally 



« 1, 



(97) 



which indicates the second term of adiabatic expansion expression Eq.f l85p can be ignored, 
or the corrections from beyond adiabatic approximations are quite small (or tiny), and hence 
the adiabatic approximation is legitimate for solving this dS-SR Dirac equation of the atom. 

Thus, we achieve an interesting consequence that the fundamental physics constants 
variate adiabatically along with cosmologic time in dS-SR quantum mechanics framework. 
As is well known that the quantum evolution in the time-dependent quantum mechanics 
has been widely accepted and studied during past several decades (see, e.g., [7]). It is 
remarkable that the time- variations of /i (or mg), h and e (see Eqs. (!77|) (!78|) ( I79|l ) belong 
to such quantum evolution effects. 
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8 2si/2-2pi/2 splitting in the dS-SR D irac equation of 
Hydrogen 

In the Section IV, we have pointed that the state of 2s^/^ and state of 2p^/^ are complete 
degenerate to all order of a in the E-SR Dirac equation of Hydrogen described in Hamiltonian 
H^ir, h, /i, e) of Eq.(j2S])- In this section, we calculate the (25*^/^ — 2p^/^)-splitting duo to dS- 
SR effects. 



8.1 Energy levels shifts of Hydrogen in dS-SR QM as perturbation 
effects of E-SR Dirac equation of atom 

The dS-SR Dirac spectrum equation for Hydrogen atom has been derived in Section VI 
(17^ . which is as follows 



where 



H^ds-SR)'^ = ( Ho{r, ht, /it, et) + H' )ij = E^, (9^ 



e2 

Ho{r, ht, fit, et) = -ihtcd ■ V + -, (99) 

r 

H' = ^{H'^ + H')^H[ + H'„ (100) 



where 



H[ = — {a^Ho{r, ht, /it, et) + Ho{r, ht, fxt, et)a^) , (101) 

Comparing above equations with Eqs. (17^ — (1761) . hereafter we have removed the subscript 
B of vb and V^, and the tilde notation ~ of for simplicity. And we have also rewritten 
the perturbation Hamiltonian H' to be explicit hermitian. In the spherical polar coordinates 
system the operator = di in Eq. fllOip is as follows 

r. d — )■_ ^ d ^ Id sin 09 

Oi = -—= I ■ V = sm6'cos0— + cos6'cos0-— ^ttt- (103) 

ox^ or r oU r sm 6 ocp 

Comparing dS-SR QM with ordinary E-SR Dirac equation of Hydrogen, there are two 
distinguishing effects in the dS-SR QM descriptions of distant Hydrogen atom (or one- 
electron atom) in the Earth-QSO reference frame: 1, The physical constants variations 
with cosmic time adiabatically, which have been discussed in the previous section (see Eqs. 
([77]) -((79])); 2, Perturbation effects arisen from H' of Eq.flTOOll in H^ds-SR) of Eq.(|98]). In 
this section, we focuss the latter. 

For adiabatic quantum system, the states are quasi-stationary in all instants. And 
hence to all instants the quasi- stationary perturbation theory works. When H(^ds~SR) = 
Ho{r,ht, fit,et) + H' , the unperturbed quasi-stationary solutions of Ho{r,ht, fit,et)ip = Wijj 
are the same as Eqs.(l33l) — (I4T1) except h ^ ht, ^ fit, e — )■ et. Then the energy levels shifts 
due to H' of Eq.f llOOp are computable in practice by the perturbation approach in QM. 
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Those shifts AE'* = Wi — E due to H' are determined by 

det {{H^ds~SR))ii' - E5u') = 0, (104) 
where i = {n, I, j, k, rrij}, H(^dS-SR) has been given in Eq. llUSi) and the elements are 

{Hids-SR))n' = mo\^) + = Wi5u' + {H'),,, (105) 

where Wi = Wn^^ are shown in Eq. (l34|) . 

Firstly, we compute the elements {H')iii = {{H[ + -?/2))m' "with i = i' . From Eq.f llOip . we 
have 



^Jdn (xrnf)a^x"^(f) - X^^i^VxT^i^)) . (106) 
Substituting Eqs.dST]) §B> into Eq. ffT06|) . we get 

1 rffixrt(f)aV:";^(f) = (107) 

I dnx":^J{r)a'x7'i^) = 0, (108) 

and hence 

{H[)ii = 0. (109) 

It is easy to be sure the validness of Eqs. fllOTp fllOSp . Considering, for instance, the case of 
state \i) = |(2pi/2)m,=i/2^ (^ggg ^ ([52])), we can calculate the left sides of flTOTl) ffTOHjl 
explicitly: 



f^^xl=J(f)(2pi/2)f^V-K=-l(f)(2pl/2) 

r f fl .Pi A/O 1 \ / A/^coserr^(e0) - ./isinee-W(e0) 

= -J-F°(^0)t,J-y,i(^0)t ( ^ M V3 I V3 1^ 

y l^V3^^'"^'V3^^'"^^VlO;\^ y^sinee^<^Ff 1(^0) + y^coseFO(^0) 
— i /"^'^ 

= — dcpsincj) d9sm9cos'^9{sm9-cos9) = 0, (110) 

y"c^fiX-i=-l(f)(2pi/2)a^Xyil(f)(2pi/2) = (^y" C?^xi=!(f)(2pi/2)C^^X-K=-l(f)(2pi/2)^ =0. (Ill) 

Therefore Egs- ffTOTj) ffTOHD hold to be true for the state of \i) = |(2pi/2)m,.=i/2^^ 
From Eq. fll02p . we vave 



(i/2). = ^(^1 (^^^c— - ^nc— j \^) = -^c{Ap V) = 0, (112) 

where = —i^, and the fact that the average value for to the stationary bound state 
\i) must be vanish have been used. Eq.f lll2l) can also be checked by explicit calculations 
based the known wave functions Eqs. fl36l) — fl38|l . 
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Combining Eq.([II2]) with Eq.dnO]), we find out that 



{H'),, = {H[)u+{H'^)u = Q, (113) 

which means {H')iii is an off-diagonal matrix in the Hilbert space. As is well known that 
the energy shifts for non-degenerate levels due to H' are as 

which could be thought as the perturbation solution of Eg. (11041) for non-degenerate case. 
Thus, since H' oc 0(l/i?2) and noting Eq.([Il3]), the non-degenerate level shifts AE'^ are 
(9(1/-R^), which are beyond the considerations of this paper. We will show in next subsection 
that the meaningful 0(l/i?^)-energy level shifts due to off-diagonal perturbation interaction 
H' are occurred for degeneration levels. Typical example is 25*^/^ — 2p^/^ splitting due to H'. 



8.2 28^/2 _ 2pi/2 splitting caused by H' 

In the Section IV, we have shown that the state of 2s^/^ and state of 2p^/^ are complete 
degenerate to all order of a in the E-SR Dirac equation of Hydrogen (see Eq. fH^ ). The de- 
generate will be broken by the effects of H'. In this subsection we calculate the 2s^/^ — 2p^/^ 
splitting caused by H'. 

By using the explicit expressions of 2s^/^- and 2j9^/^ wave functions Eqs.f l43p — ( 152|) . all 
matrix elements of H' = H[ -\- H'2 between them can be calculated, i.e, 

)2Li/2,2L'V2 = (-^l)2Li/2,2L'i/2 + (-^2)2Li/2,2L'i/2 ' (H^) 

where 

= J drr' J dn V^5i;=i/2(r) H', ij'$,,^^^,/,{r), (116) 

where {L,L'} = {s,p}, i = 1,2, and V'(2l)j=i/2('') given in Eqs.( l43|) — ( |52|) . The matrix 
element calculations are presented in step by step in Appendix C, and the results are hsted 
as follows: 

1. if[-matrix elements: H[ is given in (llOip . i.e., 

H[ = — {a^Hoir, h, /i, e) + HQ{r, H, fi, e)a^) , (117) 

where the subscript t of h, fi, e has been removed, because there is already a factor of 
(l/i?^) in the if[-expression and l/i?^-terms are ignorable. By means of straightfor- 
ward calculations (see Appendix C), we have 

,rr/vl/2,-l/2 _ ,rr/\-l/2, 1/2 _ _ / rr/ \ 1/2, -1/2 _ _ / rr' \ ^ ^2, 1/2 
\-"l/2sV2,2pi/2 — \-"l/2si/2,2pi/2 — \-"l/2pi/2,2si/2 ~ l/2pi/2,2sV2 

and others = 0, 
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where W = W(^n=2,K=±i) (see Eq. llMI) ). and notations of kc, Wc, A, s, k have been given 
in Eg. (141 1) . The above results can be expressed in exphcit matrix form as follows 



(2L^L'^'y-S- : (2si/2)V2 (2si/2)-i/2 (2pV2)i/2 (2pi/2)-i/2 

(251/2)1/2/ 

21 n n —VP), n 

(119) 



- j^2pl/2)l/2 



(2pi/2)- 



■1/2 



v 




















-iQi 




























where the matrix row's and column's indices have been labeled explicitly. Compactly, 
Eq. (lll9p can be written as follows 



Q^Q^W I kl-W^ fWc kc 



2. H2 matrix elements: H2 is shown in Eq. fll02p : 



1 "9 /i,^ „,l(9sin09 ,, 

where 7— = smt^cos07- — h cosb' cos0-— ttttt- (123) 

ox^ or r ou rsinf 



By means of straightforward calculations (see Appendix C), we obtain all elements of 



,t7-/\l/2, -1/2 _ ,„-/\-l/2, 1/2 _ _/o-/\l/2, -1/2 _ _/rr/\-l/2, 1/2 
\-"2/ocl/2 0„l/2 — \-n2/9,l/2 Or,l/2 " \-" 2/ 0„l/2 9,1/2 " \-n2/o„l/2 



'2si/2,2pi/2 \-'-'2/2si/2,2pl/2 \-'-'2/2pi/2,2si/2 V"'"' 2/ 2pi/2 ,2^1/2 

= , — ^ - 2(- + l)kc -Wc + - — 

and others = 0. 



The matrix form of (11241) is as follows 



where 



{(^2).}H,e^ ~T )' (^^^) 



02 = f# - 2(- + l)fcc - W^c + ^W-c ] ■ (126) 



2i?2 Q^AW^ - kl \Wc 's ' " kcs 

3. Since H' = H[ + H'^ (see Eq.flTOOl)). the matrix form of H' is: 

= mU'} + = ( ) , (127) 
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where 

e = Bi + 82 



4i?2 3 



ml - kl 



w, 



c 



2Wc 



> + l) 



c-Wc + ^Wl)(,128 



Obviously, e oc 0{1/R^). 



4. Substituting Eq. fll27p into Eq. fllOSp and Eq.f ll04p . we get the secular equation for eigen 
value E in the degenerate perturbation calculations: 



W -E 








-iO 





W-E 


-iO 








iO 


W-E 













W-E 



0. 



The real energy solutions are 

W-E = ±Q, orEW = H^ + e, E^-^ = W-Q, 



(129) 



(130) 



and hence 



(Ai?)(2,i/2_v/2)^i?W-i5(-) = 

'Wc 
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W 




Wl ( 






jAWl 






hcX 






Q^AWl 




Wc 



kc 

4 



2Wc 



> + 



Wc + 



kcs 



W, 



c 



(131) 



which is 0{1/R^) and the desired expression of energy level splitting of (2s^/^ — 2p^/^) 
due to H' . Eq. fll3ip represents an important effect of dS-SR to one-electron atom, and 
it is the main result of this paper. 

5. The eigenstates with eigenvalues E^^'^: Generally, the eigenstates of H' are 

= &^\2s^'Y^ + d^\2s'"^)-^'^ + d^\2p"Y'^ + &t\2p"Y^/\ (132) 
where {|2si/2)™j, |2s^/2^™^} G |r^LJ■)™^ and {i = 1,2,3,4) satisfy the following 



eigen-equation corresponding to Eq. fll29p : 

/wo -te\ f ci^^ \ 



\ 



W -iQ 
zQ W 
iB W 



J 



(133) 



Substituting f ll30p into fll33p . we get the eigenstates with eigenvalues E^^'^: 

|^(+)) = l(|2sl/2)l/2 + |2sl/2)-l/2 + ,|2pV2)l/2 + ^ 1 3^1/2^ -1/2) ^ 
\E^~)) = l(|2sl/2)l/2 + |2sl/2)-l/2 _ ,|2pV2)l/2 _ , 1 2^1/2) -1/2) ^ 



(134) 
(135) 



which satisfy = = 1, and (E(-)|E(+)) = 0. 
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6. Numerical discussions: In order to getting comprehensive indications to Eg. (11311) . we 
now discuss (A_E)(2si/2_2pi/2) = -E*-"^^ — E^~'> = 26 numerically. To Hydrogen's 2s^/^- 
and 2p^/2_g^^tgg^ ^ nie = 510998.910 eV/c^, Z = I, n = 2, k = ±1. Substituting 
them into Eq.dHl]), we get W = 510995.51 eV/c^. And then, by Eq.(!4l]), we further 
have kc, Wc, A and s. Inserting all of them into Eq. fll3ip . we obtain 



AE{z) = (AE)(2,i/2_2pi/2) 



X 358.826 eV 

X 8.36 X 10^ {Lamb shift), (136) 



where Q°(z) = ct{z) and Q^{z) have been given in Eqs.f l20|) and flMl) respectively 
(see also FIG. 2 and FIG. 3). In the expression of function AE{z) of Eq. (ll36p . when 
z were fixed, the only unknown number is R which is the universal parameter of dS- 
SR. Therefore, it is expected that R could be determined through accurate enough 
observations of the level spectrum shifts of atoms on distant galaxy. The curves of 
AE{z) of Eq. (fT36|) with \R\ = {0.5 x lO^Gly, lO^Gly, 2 x lO^Gly} are shown in 
FIG. 6. It is essential that FIG. 6 shows that when z > 1, AE{z) » 1 Lamb shift. 
This indicates that comparing with the usual QED's hyperfine structure effects (i.e., 
the Lamb shift measured in the laboratory), the dS-SR QM fine structure effects are 
dominating for the splitting between 2s^/^- and 2p^/^- states of Hydrogen atom on 
distant galaxy. In the TABLE I, the AE{z) for \R\ = {lO^Gly, lO^Gly, lO^Gly} and 
2; = {1, 2} is listed. It is learned that AE{z) » 1 Lamb shift to all cases too. 



2 



C3 




Figure 6: Functions AE{z) of Eq.dHI]) with \R\ ^ 0.5 x IQ^Gly, lO^Gly, 2 x lO^Gly arc shown. The unit 
of the energy spUtting AE{z) is (Lamb shift)~ 4.3 x 10~^eV. 



9 Summary and discussions 

In this paper, we have solved the de Sitter special relativistic Dirac equation of Hydrogen 
in the Earth-QSO framework reference by means of the adiabatic approach and the quasi- 
stationary perturbation calculations of QM. Hydrogen atoms are located on the light cone of 
the Universe. FRW metric and ACDM cosmological model are used to discuss this issue. To 
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Tabic 1: The energy level splitting of Hydrogen's {2s^^'^ - 2p^/'^), AE{z) = (Ai;)(2si/2_2pi/2) (see Eqs. 



(|13ip (|136p ): R is the radius of de Sitter pseudo-sphere in dS-SR, which is a universal parameter in the 
theory and \R\ > Rh = 13.7G^?/(Horizon of the Universe). Gly = 10^ light years, z is the red shift. (Lamb 
shift)- 4.3 X IQ-^eV. 



\R\ 


W^Gly 




lO^Gly 


z 


1 2 


1 2 


1 2 


AE{z) (eV) 
(Lamb shift) 


2.6 11 
6 X 10^ 2.7 X 10^ 


2.6 X 10^2 0.11 
5977 26541 


2.6 X 10^^ 1.1 X 10"^ 
59.77 265.41 



the atom, effects of de Sitter space-time geometry described by Beltrami metric are taken 
into account. The dS-SR Dirac equation of Hydrogen turns out to be a time dependent 
quantum Hamiltonian system. We have provided an exphcit calculation to examine whether 
the adiabatic approach to deal with this time-dependent system is eligible. Since the radius 
of de Sitter sphere \R\ is cosmologically large, it makes the time-evolution of the system is 
so slow that the adiabatic approximation legitimately works with high accuracy. Based the 
dS-SR Dirac equation's solutions up to 0{1/R^), some remarkable effects of dS-SR atom 
physics are revealed: 

1. The fundamental physics constants variate adiabatically along with cosmologic time in 
dS-SR quantum mechanics framework. As is well known that the quantum evolution in 
the time- dependent quantum mechanics has been widely accepted and studied during 
past several decades. It is remarkable that the time- variations of fi {or rrie), h and e 
(see Eqs. ( 1771) ( 1781) ( 1791) ) belong to such quantum evolution effects. 

2. The fine-structure constant a = e^/ (he) keeps invariant along with time up to (9(l/i?^) 

In the expression of a, the e^'s time- variation and Ws are canceled each other. 
However, whether or not this cancelation mechanism works to the next order 0{1/R'^) 
remains to be open so far. 

3. (2s^/^ — 2p^/^)-splitting due to dS-SR Dirac QM effects: Distinguishing from E-SR 
Dirac QM theory of Hydrogen atom, the degeneracy of (2s^/^ — 2p^^'^) is broken in 
dS-SR QM. By means of the quasi- stationary perturbation theory, the (2s^/^ — 2j»^/^)- 
splitting AE{z) has been calculated analytically, which belongs to (9 (1/i?^) -physics 
of dS-SR QM. Numerically, we found that when \R\ ~ {lO^Gly, WGly, lO^Gly } 
(note the Universal horizon Rfj ~ 13. 7 Gly « \R\), and z ~ {1, or 2}, we have 
AE[z) » l(Lamb shift). This indicate that for this case the hyperfine structure 
effects due to QED could be ignored, and the dS-SR fine structure effects are dominant. 
Therefore, we suggest that this effect could be used to determine the universal constant 
R in dS-SR, and be thought as a test of new physics beyond E-SR. 

Finally, we address again that the dS-SR is a natural extension of E-SR. What we achieved 
in this paper are that we revealed new effects in one-electron atom physics, which are beyond 
E-SR and hence can be used to recognize dS-SR. 
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Appendix A: Electric Coulomb Law in QSO-Light-Cone 
Space 

Let's derive (!62|) . The action for deriving electrostatic potential of proton located at = 
{Q° = ct, , = 0, = 0} within background space-time metric g^^ = B^^{Q) of 
eq.( fT6l) in the Gaussian system of units reads 

S = - I d'L^F.^F^'' --J d'L^gfA^, (137) 

where g = det{B^^), F^iy = ^ - ^ and = {f = cpproton/ VB^, j} is 4-current 
density vector of proton (see, e.g, Ref. [27] : Chapter 4; Chapter 10, Eq.(90.3)). The explicit 
matrix-expressions for B^y{Q) and B^'^{Q) up to 0{1/R^) are follows: 



/ 1 



2(Q°)^-(Qi) 

i?2 



/?2 





Vo 












(Q0)2-(Ql)2 








ij2 



Bi^^iQ) 



i?2 



-1 - 





vo 



2(qV-(Q»)^ 






-1 + 







/?2 



o\ 










— w 



Making space-time variable change of L^^ ^ — = = = ct^ — ct, = L* — Q*}, 
we have action S as 



S 



IGvrc 



- I d^x3^{x)A^{x)\ y/-det(5^,(Q)) 



:i38) 



and the equation of motion 5S/5A^{x) = as follows (see, e.g, [27], Eq. (90.6), pp257) 

in 



d^F^"" = B^^d^F^^ = -—f. 



(139) 



In Beltrami space, A^^ = {0b, A} (see, e.g., [27], eq.(16.2) in pp. 45) and 4-charge current 
= {cpproton / VBqo, j}. Accordiug to the expression of charge density in curved space in 
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Ref. [27], {pp.256, Eq. (90.4)), Pproton = Pb = ^S^-^K^) and j = 0, where 

7 = det(7ij), 



dp = 'jijdx'^dx-' 



Noting Bqi = Bio 



BpiBjo 



00 



goigjo \ 

9oo ) 

dj'JC diijC 



dx'dx^ {see eg. (84.7) in Ref. [37]) 



, Boo ~ 1, and fioi^io - 0{1/R^) ~ 0, we have 



and hence 



Vl = \/det(7ij) - ^-det{Bij) 
e5(3)(x) 



Pproton — PB 



3 = 0. 



^/^det{B~) 

1. When = in Eq. (11391) . we have the Coulomb's law (1601) . i.e., 



(140) 
(141) 

(142) 
(143) 



B'\Q)did,(PB{^) 



{Qr - {Q 
— 47re 



1\2 



{Q 



1\2 



i?2 d{x^Y 



^b{x) = 



= 5(3) (. 



(144) 



^-det{B,^{Q))Boo{Q) 

where B^Q) = + ^^^^^^.i^.i + ^^^^5..^,. + ^^^^^^.s^.s + 0{R-^ 
has been used, and Bij were given in ( !T6|) . Expanding (11441) . we have 



+ 



92 ^2 

+ 



X ^1- 

= — 47re 



2(g")^ - { Q 

2i?2 



1\2 



1 + 



An 1 



3{Qr-A{Q 



1\2 



{Qy - {Q 

R^ 



1\2 



2i?2 



e6{x^)6{x^)6{x^ 



1 - 



Noting 5(x^) = (5(xV[l + {Qy/2R'^]){1 - {Qy/2R^), we rewrite above equation as 
follows 



^2 



92 92 
+ 77^^ + 



+ (g)!])2 a(a;2)2 9(x3)2 



fx) = -Aire 1 



my - i Q 
2i?2 



l\2l 



X 1 



^^'^'^5(a:V[l + ^])5(xW) 



2i?2 



-Aire 1 



3(Q0)2 _ 2{Q 



1\2 



2i?2 



5(a:V[l 



(Q 



1\2 



2i?2 



-])6{xy{x' 



Setting 



{Q 



1\2 



the above equation becomes 

^2 



92 92 
+ 



9(xl)2 9(x2)2 9(x3) 



-47re 1 



3(g°)2-2 (g 
2i?2 



1\2 



(145) 

6{x^)6{x^)6{x'). 

(146) 
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Then the solution is (pnix) = (l - ^Mli^HMlij e/rs with 



Therefore, we have 

3(<3Y-2(Q')^\ e 



which is the scalar potential in Eq. (l62l) in the text. 
2. When /i = i (i = 1, 2, 3) in Eq. ffT39|) . we have 

d'df.A^' - B'^'d^d^A' = -^f = 0. (148) 

By means of the gauge condition 

d^A>' = 0, (149) 

we have 

B^'^d^d^A' = 0. (150) 

Then 

A' = (151) 

is a solution that satisfies the gauge condition fll49p (noting OqA'^ = -^4>b{^b) = 
due to 1^ = If^ = ). Eq. (fT5B is the vector potential in Eq.([62]) in the text. 

Appendix B: Adiabatic approximative wave functions in 
<5^ci?-Dii*ac equation of hydrogen 

Now we derive the wave function of (185!) in the text. We start with eq.( l82|) . i.e. 

ihdtiJ = H{t)ij = [Ho{rB, h, fx, e) + //^(t)]^^, (152) 



where 



H{t) = Ho{rB,h,fi,e) + H'o{t), (153) 



e2 



Ho{rB,h, fi,e) = —ihca-VB + l^'C^P (gee eg. (l33|) ) (154) 

tb 

m) = -(^^^Ho{rB,h,3fi,e). (155) 

Suppose the modification of H(t) along with the time change is sufficiently slow, the 
system could be quasi- stationary in any instant 6. Then, in the Shrodinger picture, the 
quasi- stationary equation of H{6) 

H{d)Un{^, 9) = W„(^)t/„(x, 9) (156) 
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can be solved. By fll53p fll54p fllSSp and t 6, the solutions are as follows (similar to eq. 
in text) 



Wn{9) = WnAd) = /i9cMl + 



a 



n — \k\ + sy 



-1/2 



(157) 



p2 2 

hec he ^ i?4 



|«:| = (j + l/2) = l, 2, 3 



s = V k'^ — a^, n = 1, 2, 3 
where (see dZZD, ([78]), ([79]) in text) 



he 



2i?2 

2i?2 
1 - -TTTT I e, 
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The complete set of commutative observable is {H, k, j2, j.,}, so that we have 

Un{yi,0) = tpn,K,j,jA^B,he,He,eg), 

where j = L + |S, hn = /3(S - L + h). [t/n(x, 6^)] is complete set and satisfies 

j £xUn{x,9)U^{x,9) = 5mn, Tl = {nr,K,j,jz}. 



(158) 
(159) 
(160) 

(161) 

(162) 



Thus, the solution of time- dependent Shrodinger equation (or Dirac equation) f ll52p can 
expanded as follows 



V'(x,t) = VC„(t)f/„(x,t)exp \-i [ ujn{e)de 
Substituting (I163P into ( ]152p . we have 



Wn{e) 
h ' 



(163) 



jfi V'(C'„f/„ + C„f7„) exp -i j 

V Jo 



ujn{d)dd 



(164) 



By multiplying [/^exp 
using fll62p . we have 



j^ujm{d)d9 to both sides of eq. (ll64p . and doing integral to x by 



Cm + Cm / d^xU^Um 



' C„ / d^xU^jjn exp -i / (a;.„ 
„ J V Jo 



Um)d6 



0, (165) 



m =1,2, 3, 



where means that n 7^ m in the summation over n. Noting f ll62p . we have 



U*^Umd'x + I U:,Umd'x = 



m'-^rn^ 



(166) 
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and hence 



is purely imaginary number. Denoting 



OLr, 



U^Und X, and COnm = OJn- ^m-, 



then eq. fllGSp becomes 

Crn + i^Cyn + J]] ' C'nttmn exp -% / 

To further simphfy it, we set 

K(x,t) = t/„(x,t)exp 



0. m =1,2,3, 







then 



V^(x,t) = 5^G„'(t)K(x,t)exp -z / 

L ^0 



UJn{e)dd 



where C^{t) = C„(t)exp i J^^ (3n{6)d9\, and 



Substituting f lT72ll into ffT69ll . we finally get 



L ^0 



de 



0. m =1,2,3, ■■■ 



where 



1, 



Now let's solve fll73p . Firstly, we derive amn- By f ll56p . we have 

dH. 



dt 



By multiplying U*^ and doing integral over x, we have 



U*HUnd^x 



U*HUnd^x = Wn / WU^d'x 



T* TT j3„ 



i.e.. Hum ~l~ ^^m^mn ^^n^mni 



SO that 



a 



— , U* U d^r 



Wu-Wm 



Therefore eq. fll73p becomes 



0. m = l,2,3. 
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Suppose in the initial time the system is in s-state, i.e., C„(0) = C„'(0) = 5ns- To adiabatic 
process, H{t) — > 0, then the 0-order approximative solution of eq. fll78p is 

[C2t)]o = 6ms- (179) 
Substituting (11791) into fll78p . we get the first order correction to the approximation 



[CL]i = ^^^^P [-^I^Lde ) = 0, niy^s. (180) 



Since the dependent on time t of f/„(t) is weak for adiabatic process, eg. (11671) indicates /3„ 
is small, and by (I174p . we have u'^^ ~ UJms- Then, from (I180p . the first order correction to 
the solution is 

= {e'--' - 1) , rnj^s. (181) 

Substituting (11801) (I18ip into (11711) and neglecting we get the wave function as follows 

^(x,t) ^ f/.(x,t)e-^ + E #^ i^'^"-'' - 1) UU^,t)ei-^fo^''). (182) 

By using eqs. (11610 . (11590 . (11580 . (11600 . we finally obtain the desired results 



^s '^rns 

183) 



where 



ht=[l- ^ j h, (184) 

e,= (l-g)e, (186) 
ffTMD ffTSS]) and ffT86|) are the equations ([TTD, ([78]) and ([79]) in the text. Eq. ffT83|) is just 



Eq.(|85D in the text. 

Appendix C: Calculations of elements of the perturba- 
tion Hamiltonian H'-matric in (2s^/^-2p^/^)-Hilbert space 

Now we derive Eq. (lll8p and Eq. (11240 . We start with the dS-SR Dirac spectrum equation of 
Hydrogen, which has been shown in Eqs. (I98l) - (ll02p in the text: 

Hids-SR)i' = ( Ho{r, ht, fit, et) + H' )ij = E^, (187) 

where 

Ho{r, ht, fit, et) = -ihtcd ■ V + fitc^(3 (188) 

r 

H' = ^{H'^ + H')^H[ + H'„ (189) 
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^1 = — Jm' ("^^ol?^, K e) + HQ{r, h, fx, e)a^) , (190) 



where 

The definition of if'-elements in the i^o-eigenstate space, (i^')2Li/2 2L'i/2' been given in 
Eqs. flllSp f lll6p . The eigen values and eigen states of Hq are given in the section IV. 
(I) iJ{-matrix elements: 



1. m 



,\l/2, -1/2 . 
2si/2,2pl/2- 



^ I ^{2pi/2)(^)X«^^^(f)(2pl/2) 

V/(2pV2)(r)x:y^(f)(2pi/2)^ 

'^^B^Jo ^^^^ |^{2sV2)(r)/(2pi/2)(rydfixK^^"^(f)(2.i/2)C^^Xli (f)(2pi/2) 

-/(2si/2)(r)^(2pi/2)(r)/rffix-K"^(f)(2.i/2)C^^XK^^^(r)(2pV2) }> (192) 



where IV = W(^n=2,K=±i), and the explicit expressions of 2s^/2_ g^y^d 2p^/2 -wave func- 
tions of {^(2si/2)('^), /(2si/2)(r), ^(2pi/2)(r), /(2pi/2)(r), ^^^^(f) (2.1/2), X±y^(f )(2pi/2)} 
are given in Eqs.f l43p — f l52|) in text. From them, we have 



d^Xn"^^ (f ) (2s^/2)(r\Jj^ (f) (2pi/2) 
= l^'\4> sine (YlO) ( J J ) 



2- ^0 l\l Jfcos^Ff ^(^0) - ->/|sin^e-^'^Fi° 



I sin ee'^Y^-\e(f)) - a/ | cos 



1 1 

c/^sin^(sin^cos^ + cos2^) = -; (193) 



2io ^ '3 

d^X^^K (f ) (2s1/2)^^Xk (f ) (2pi/2) 

1 \ / V 3-^1 



^ d(j) de sin (- cos ^^0°. - sin ^6-^*^0°) ( J 



3 

l\^0 



-- [ desinOcos'^e = --, (194) 
2 Jo 3 



where Fi°(^0) = y^J^cos^, y/(e(/') = -y'|;e*<^ sin ^, Y^^^dcf)) = y'l^e"^'^ sin ^ and 
^^0(^0) = .fT have been used. Substituting Eqs.([l93]) ffTMl) into f fT92|) . we get 



(^()Jf?;2,2pV2 = Jo ^'"^^ {^(2sV2)(r)/(2pi/2)(r) + /(2si/2)(r)^(2pi/2)(r)} (195) 
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Inserting the explicit expressions of radial wave functions fi'(2si/2)('"), f{2s^/^){r), and 
5'(2pi/2)(r), /(2pi/2)(r) (i.e. Eqs.dH]), (iS]), dH]), ([50]) in text) into the integral in 
Eg. (1 195 p . and accomplishing the calculations, we have 

drr^ {£/(2si/2)(r)/(2pi/2)(r) + /(2si/2)('^)^(2pV2)(r)} 

where formula Jp°° (irr*^"^ exp(— /ir) = r(z/)//i'' were used. Consequently, substituting 
Eq. ffT96|) into Eq. ffT95|) . we obtain 



\-"l/2si/2,2pi/2 — '^^1 

which is just desired result of Eq. OllSp . and all above calculations have been checked 
by the Mathematica. 



2. By means of similar calculations we get also that 

> -1/2, 1/2 
'2sl/2,2pi/2 



{H[)-l!^^lil, = -^el. (198) 



Since H[ = H'^ , we have 



iH[)-Jiiil. = mfS)^''lr = ^el, (199) 



'2pl/2,2si/2 — V\-'-'l/2sl/2,2pl/2. 

= mrjfl\2il^r = ^©i- (200) 

Furthermore, to all other elements of iJ(-matrix, since J^^ d(j) exp(±m0) = and 
Jq dO sin cos^""*"^ = and etc, the explicit calculations show that all those H[- 
matrix elements vanish. Consequently, all elements of H[ are calculated, and Eq.f lllSp 
is proved. 

(II) ifg-matrix elements: 

if 2 hcis been given in Eqs. f ll22p fll23p in the text, which is as follows 



d 1.0 sin (j) 

where 7-— = (9i = sin6' cos0— — h cos6' cos0-— ttttt- (202) 

ox^ or r oO r sm ocp 



We derive Eq. (11240 in text now. 
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1 /if' \ 1/2, -1/2 . 



,rrnl/2,-l/2 _ 2 f JjV^^ {v) H' ^'^'"^ 

\-"2/2sl/2,2pl/2 - / "-^^ / 'A'(2s)j=l/2'^'^>' -"2 '/^(2p)j: 



(2p)i=l/2l''J 



4i?2 



nc y"cirr^o?^^|c/(2si/2)(r)[<9i^(2pi/2)(r)]x«/^Hf)(2si/2)X-y^(f)(2pi/2) 



^(2../2)(r)x;:;/^^(f)(v/^) 

V/(2pi/2)(r)x-y^(f)(2pi/2)^ 



^{2sV2)(r)^(2pi/2))(r)x«/^nr){2sV2)[<9ixiy^(f)(2pi/2)] + /(2sV2)(r)[9i/(2pi/2)(r)]x«/^nr)(2si/2)Xiy^(f)(2pi/2) 

/(2si/2)(^)/{2pi/2))(^)X?c^^"^(f)(2si/2)[^lXiy^(f)(2pl/2)] - [5l^(2si/2)(^)]^(2pi/2))(^)X«^^"^(f)(2sl/2)Xiy^(f)(2pi/2) 
^(2si/2)(^)^(2pi/2))(^)xy^"^(f)(2si/2)[9lX:y^(f)(2pi/2)] - [5l/(2.i/2)(^)]/(2pV2)(r)xy^"^(f)(2,i/2)Xiy^(f)(2pi/2) 

/{2si/2)(r)/(2pi/2)(r)xy^"^(f)(2si/2)[<9iXiy^(f)(2pi/2) 



4/?2 



5^ ^?(2pV2)(r) (?(2,i/2)(r) + /(2pi/2)(r) 



■9/(2^1/2 )(r) 

Or 



/(2si/2)( 



r 



rf^^ [^(2si/2)('^)^(2pl/2)(r) - /(2si/2)('^)/(2pi/2)(^)] f ,(203) 



where the integrals to dVL have been accomphshed in terms of the exphcit expressions 
of x±f (f)(2si/2), xiy^(f){2pi/2) in Eqs.dM]) gT]) (Ell) ([52]) and Eq.(l202]). Substituting 
expressions (jSj) ( H5l) ( H9l) ( l50l) into Eq. (12031) . and finishing the integrals, we get 

/rr/\l/2, -1/2 _ .p, 
\^2/2si/2,2pi/2 = -^'-'2 

2fl2 6v'4H'^ - itj VlVc kcs ^J' ^ ' 

which is just Eq. (ll24p . and all above result have been checked by the Mathematica 
calculations. 



2. By means of similar calculations we get also that 

>-l/2, 1/2 
'2si/2,2pi/2 



(^2)oy/2'o^/2 = -^62. (205) 



Since H2 = H2 , we have 



(^2)2„l/2'252 = ((^2)2S2"2/'/2)* = ^62, (206) 



'2pl/2,2si/2 — U-'-'2/2sl/2,2pl/2j 

1/2, -1/2 _ ,,r.,v-l/2, 1/2 ^ 
'2pl/2,2si/2 ~ V\-"2/2sl/2,2pl/2j 

-277 



iH2)l%-'!l. = m);!il'i!'.r = ^q^. (207) 



Furthermore, to all other elements of if2"niatrix, since J^^ dcp exp(±m(/)) = and 
Jq dO sin 6 cos^""'"'^ ^ = and etc, the explicit calculations show that all those H'2- 
matrix elements vanish. Consequently, all elements of H2 are calculated, and Eq.( ll24p 
is proved. 
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